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The Falicov-Kimball model with a correlated-hopping interaction is solved using an extended dy- 
namical mean-field theory that becomes exact in the limit of large dimensions. The effect of corre- 
lated hopping is to introduce nonlocal self-energy components that retain full dynamics as D — » oo, 
thus introducing an explicit fc-dependence to the single-particle self-energy. An explicit solution for 
the homogeneous phase at D — 2 reveals significant nonlocal dynamical contributions in the physi- 
cally relevant regime of a moderately large correlated-hopping amplitude, indicating that important 
nonlocal correlations are omitted in Hubbard-like models upon neglecting the correlated-hopping 
interaction. 
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The development of the limit of infinite dimensions j3 or 
the Dynamical Mean-Field Theory (DMFT), has opened 
exciting new possibihties in the study of correlated 
electron systems.0 One notable example is the Mott- 
Hubbard metal- insulator transition realized in V2-j,03 
and Cai_xSra;V03, where the main features of the opti- 
cal conductivity and the photoemission data are success- 
fully accounted for within the. DMFT treatment of the 
single-band Hubbard modcLHu This approach, however, 
fails to explain the moderate mass enhancement and the 
apparent fc-dependence of the single-particle, self-energy 
close to the transition in these materials,Q'l3 indicating 
that important nonlocal interactions are either omitted 
in the Hubbard model, or mistreated in the DMFT. 

One nonlocal interaction term clearly absent in the 
Hubbard model is that of correlated hopping (CH), in 
which hopping between neighboring lattice sites depends 
on the occupancy of the opposite spin orientation. Such 
an interaction term is an integral part of the Coulomb 
repulsion between electrons, and with estimates rang- 
ing from about 0.5eV in transition metalsjj to 0.8eV in 
the cuperatesjj and up to 3.3eV in benzenejj it is cer- 
tainly comparable to, if not larger than, the correspond- 
ing amplitude for single-particle hopping in real systems. 
Despite this fact and in spite of Hirsch'|S-|Suggestion of 
a new mechanism for superconductivity,t3 there is still 
no qualitative understanding of the effects that a CH 
term might have. While some rigorous statements can 
be made, these aije-|ijeatricted to certain models in one 
spatial dimension,EirM ^d to a few very special cases 
in higher dimensions .EjEEI There is thus an obvious need 
for a careful investigation of CH within a reliable non- 
perturbative approach, applicable to the entire range of 
interactions from weak to strong coupling. 

In this paper, we generalize the limit of infinite di- 
mensions, D oo, as to include the jSH interaction. 
Focusing on the Falicov-Kimball modelO — a simpli- 
fied Hubbard model with one spin species tied down 
— the effect of CH is to introduce nonlocal self-energy 
components that retain full dynamics as D ^ oo. This 
marks departure from the local approximation. Similar 



to existing formulations of the infinite-dimensional limit, 
though, both the local and nonlocal self-energy compo- 
nents are extracted from an effective single-site action, 
which can still be treated analytically for the Falicov- 
Kimball model. Solving the resulting DMFT for the ho- 
mogeneous phase at D = 2 and in the physically relevant 
regime of a moderately large CH amplitude, the dynam- 
ical contribution of the nearest-neighbor self-energy is 
found to be comparable to that of the local self-energy, in- 
dicating that important nonlocal correlations are omitted 
in Hubbard-like models upon neglecting the CH interac- 
tion. This not only suggests that CH may be essential for 
the complete understanding of the metal-insulator tran- 
sition in V2-y03 and Cai-ajSr^jVOa, but further calls 
into question its omission in any Hubbard-like modeling 
of real systems. 

We begin our discussion with the Hamiltonian 

<ij> 

+ {Ef-^,)Y,f^f^'~^^Y.'^\d^, (1) 

i i 

in which itinerant d fermions interact via an on-site 
Coulomb repulsion U with localized / fermions on a D- 
dimensional hypercubic lattice. Motion of the d fermions 
has two components: a single-particle, nearest-neighbor 
hopping term ti, and a correlated-hopping term t2, in 
which the hopping amplitude for the d fermions depends 
on the occupancy of the corresponding / fermions. The 
two fermion species share a joint chemical potential, /i. 

For la = 0, Eq. (|l|) reduces to the Falicov-Kimball 
model,EZI which is perhaps the simplest many-body 
Hamiltonian to exhibit both a Mott-Hubbard metal- 
insulator transition and long-range ordered phases. A 
nonzero t2 modifies the single-particle dynamics of the 
itinerant fermions in two ways. First, there is a straight- 
forward Hartree renormalization of the static single- 
particle hopping amplitude according to ti ^ t = ti + 
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FIG. 1. Representative d self-energy diagrams. Here, 
squares indicate a correlated-hopping vertex (t2); full 
(dashed) lines represent bare d (/) propagators; and j and 
k are distinct nearest neighbors of site i. In Fig. (a), there is 
summation over all equivalent sites k. The diagram of Fig. (a) 
gives a contribution of order D~^^^ to Eij, while the diagram 
of Fig. (b) gives a contribution of order to Efcj. Both 

terms are relevant in the limit D ^ oo. 



2n ft2 , where Uf is the / occupation number (assuming a 
homogeneous phase). In addition, the d and / fermions 
acquire new self-energy terms with full dynamics, which 
survive also in the absence of an on-site U. To investigat 
the effect of CH we employ the limit of large dimensions,! 
which has proven useful in studying the Falicov-Kimball 
model.ESo To this end, both the single-particle and CH 
terms in Eq. (|^) have been properly rescaled with D. 

For conventional Hubbard-like models, the self-energy 
is purely local at D = oo, rendering the local approxima- 
tion exact. This paradigm breaks down as soon as CH is 
considered, as exemplified in Fig. 1 for the Hamiltonian 
of Eq. (0). In Fig. 1, sites j and k are distinct nearest 
neighbors of site i. Thus, after summation over all equiv- 
alent fc, the diagram of Fig. 1(a) gives a contribution of 
order to the nearest- neighbor d self-energy I]y(z), 

while the diagram of Fig. 1(b) (in which both j and k 
are fixed) gives a contribution of order 1/D to the next- 
nearest-neighbor self-energy Sfcj(z). Since both orders in 
D coincide with the necessary large-D rescaling of the 
corresponding nearest- and next-nearest-neighbor hop- 
ping amplitudes both self-energies contribute in the 
limit D ^ oo. As shown below, the nearest- neighbor and 
next-nearest-neighbor self-energies are in fact the only 
nonlocal self-energy components relevant at Z? = oo. 

A systematic formulation of the D oo limit for the 
Hamiltonian of Eq. (^ employs the Luttinger-Ward func- 
tional $ — the sum of all vacuum-to-vacuum skeleton 
diagrams. Consider a generic $ diagram. Labeling each 
vertex according to the site index of its incoming and 
outgoing / lines, we attach the I/a/D factor associated 
with each CH vertex to its incoming or outgoing nearest- 
neighbor d line. With this convention we notice that, 
even after summation over all intermediate sites, any 
path connecting two vertices at sites i ^ j gives a con- 
tribution of order D""'^-'!!/^ or higher to the diagram at 
hand (||i — j|| is the minimal number of nearest-neighbor 
steps leading from site i to site j). Given that (i) any two 
vertices within a $ diagram can be connected by four in- 
dependent paths having no lines in common, and (ii) the 
total number of equivalent sites j scales as dH'"-?!!, the 
overall contribution of such a diagram (after summation 



over all equivalent j's) is O (Z3~"*~^ll). Thus, to order 
Z3~", only skeleton diagrams with a maximal inter- vertex 
distance \\i ~ j\\max1^ n need to be considered. Specifi- 
cally, only purely local vertices are left at D ~ oo. 

While the above classification appears identical to that 
for conventional Hubbard-like models, it should be em- 
phasized that purely local vertices in our convention do 
not imply purely local diagrams. To the contrary, a CH 
vertex at site i couples the operators di and fi to the 
nearest-neighbor shell operator 
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in which j is summed over the 2D nearest neighbors of 
site i. Thus, a CH vertex at site i necessarily introduces 
at least one propagator that is not confined to site i. 

Despite the above nonlocality, the D oo limit re- 
mains tractable since $ decouples in this limit into 



(3) 



where $;oc is the functional generated by the action 
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Here Q{t — t') is a 2 x 2 matrix propagator, chosen such 
that, when dressed with the effective-action self-energy, 
it coincides with the corresponding lattice propagator: 



with 



G{iLo„ 



Gr(*c.„), Gr(^c.„) 



(5) 
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[ujn = 7r(2n + 1)T are the Matsubara frequencies]. Here, 
in using the 2x2 matrix notation of Eqs. (||)-(|6|), we 
have identified the indices a = dyip with a — 1,2, re- 
spectively. Contrary to 5, which requires knowledge of 
the fully dressed lattice d propagator, the bare / prop- 
agator in Eq. (^ is simply Qj^{iujn) — i^n + M ~ ^/j 
which stems from the localized nature of the / fermions. 

From Eq. (||) it is apparent that $ at Z? = oo is a 
functional of the local d and / propagators, as well as 
the nearest-neighbor and next-nearest-neighbor d Green 



functions which enter via G^fjCff", and Gj^^. Thus, in 
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addition to the local d and / self-energies, one must also 
consider the nearest- and next-nearest-neighbor d self- 
energies. From a functional derivative of Eq. (||) one has 
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(7a) 
(7b) 

(7c) 



where «i,j::>> denotes distinct lattice sites with a com- 
mon nearest neighbor. In Eq. (|7c|) , -/V«ij» is equal to 
one for i and j on the same axis; otherwise it is equal to 
two. Accordingly, the fc-dependent d self-energy reads 



Sg(z) = I]f,(z) -f eg [I]f^(z) + S^,(z)] 



(8) 



where eg 



\/2/D cos(fc;). Notice that, contrary 
to conventional formulations of the limit of large dimen- 
sions. Eg has an explicit k dependence which enters solely 
through the dimensionless energy eg: ^^{z) = 'S^giz). 

A self-consistency loop can finally be closed by express- 
ing the lattice Green function of Eq. (||) in terms of Sg. 
For the spatially homogeneous phase this gives 



G(z) 



de 



z + /i — iie — Ef (z) 



(9) 



where p(e) is the density of states for the energy eg, i.e., 
a Gaussian for D = oo. Equation (^) can further be writ- 
ten in closed form in terms of the bare (U,t2 = 0) tight- 
binding d propagator. We also note that Eqs. (||)"(H) 
were derived under the explicit assumption of a hypercu- 
bic lattice with only nearest-neighbor single-particle hop- 
ping. Different self-consistency equations will generally 
apply to other lattices and other tight-binding models. 

Equations (^)-(^ are an exact formulation of the limit 
D ^ oo for the Hamiltonian of Eq. (|l|). In particular, 
for t2 — 0, when only S^^(ict;„) is nonzero, they properly 
reduce to the corresponding equations for the Falicov- 
KimbaU model at D = c».ll3 In the spirit of the DMFT, 
the same set of equations can also be applied to any finite 
dimension provided the actual i'-dimensional density 
of states p{e) is used in Eq. (^). Below we present results 
for D = 2. Although designed at present to describe the 
homogeneous phase, this formulation is easily extended 
to phases with long-range order, as well as to other lattice 
models with CH. For example, in the Hubbard model 
one has to introduce separate tpia- fields for each spin 
orientation, but the basic formulation remains the same. 

A crucial advantage of the Hamiltonian of Eq. (|l|) over 
its Hubbard counterpart is that the action of Eq. H) can 
be treated analytically. Integratiiig out the / particles, 
the lattice Green function of Eq. m) is expressed as 
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FIG. 2. The d spectral function as a function of t2/ti, for 
Uf + Ud = 1, U/ti = 4, T/ti = 1, and D = 2. For both 
t2/t\ = Q and t2/t\ = —1, the half- filled model is particle-hole 
symmetric, with a Mott gap for U /ti — 4. This gap gradually 
closes upon increasing t2/ti > or decreasing t2/ti < —1, 
until the lower and upper Hubbard bands are merged into 
one broad band. Note the general resemblance between curves 
with values of t2Ai that are symmetric about t2/ti = —1/2. 



G{iuJn) = {l-nf)g{iu!n) + 71/ 
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where n/ is the / occupation number: rif = l/{e"'^ + 1), 
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Ef ~ + - ^ Indet 
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U t2 
t2 
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Once Eqs. (||), (P)-(|lT|) have been iterated and the / 
occupation number determined, the same equations can 
be used to calculate G{z) for an arbitrary complex z. 

Restricting attention to D = 2 and Ef = 0, we fo- 
cus hereafter on the homogeneous phase of the half-filled 
case: +rtrf = 1, where Ud is the d occupation number. 
Setting U/h = 4 and T/ti = 1, Fig. 2 shows the evo- 
lution of the d spectral function as a function of ^2/^1- 
For the Falicov-Kimball model, t2 = 0, the d spectrum 
has a Mott gap. This gap persists throughout the range 
— 1.4 < t2/ti ^ 1, although the chemical potential does 
not always fall inside the gap. For either t2/ti > 1 or 



t2/t: 



< 



-1.4, the lower and upper Hubbard bands merge 



into one broad band, reflecting the increase in kinetic en- 
ergy of the itinerant fermions. Due to the lack of particle- 
hole symmetry, the d spectrum is generally asymmetric. 
Indeed, only for t2 — and ^2/^1 = —1 is the half- filled 
model particle-hole symmetric, in which case the result- 
ing spectrum is temperature independent .Ell 

As shown in Fig. 3, one cannot simply absorb t2 
into an effective Hartree renormalization of the single- 
particle hopping within the simpler Falicov-Kimball 
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FIG. 3. The d spectrum for the Hamiltonian of Eq. (|l|) at 
D = 2 versus the one obtained within the DMFT for the 
Falicov-Kimball model with the renormalized single-particle 
hopping, tpKM ~ ti + 2t2nf. Here n/ and Ud = 1 — n/ are 
kept the same in both models, as are U jtx and T/ti. 



model. Specifically, tuning /i and Ef in the latter model 
as to maintain the same and n^^^ there are substantial 
deviations between the d spectral function for the Hamil- 
tonian of Eq. (|^) and that obtained within the DMFT for 
the Falicov-Kimball model with the renormalized hop- 
ping, tpxM = ^1 + 2t2'T'/- The d spectrum for a nonzero 
^2 is notably more asymmetric, with considerable broad- 
ening (narrowing) of the upper (lower) Hubbard band. 

The effect of CH is best seen, though, in the nonlo- 
cal d self-energy components, which are absent in the 
DMFT for t2 = 0. Figure 4 depicts the imaginary parts 
of Eo, Si, and E2, defined from the expansion Sr: = 



defined from the expansion 
So + Siej + S2e2. From Eqs. {^~§), Si = %/2^S<,j> 
and E2 — 2£)I]«ij»/A^«i,j» are the appropriately 
scaled nearest-neighbor and next-nearest-neighbor self- 
energy components. Evidently, the imaginary part of 
El (which has no Hartree contribution) is comparable to 
that of Eq in the physically relevant regime, 1^2/^1! ^ 1- 
Hence, important nonlocal dynamical contributions are 
omitted in the DMFT upon neglecting ^2- 

Anticipating a qualitatively similar result for the Hub- 
bard model, CH thus provides a natural mechanism for a 
strong /c-dependence in the single-particle self-energy, in- 
cluding for D — > CX3. This may prove important, e.g., for 
the understanding of the apparent fc-dependence and the 
moderate mass enhancement close to thejuetal-insulator 
transition in V2-j,03 and Cai-ajSra^VOaJj^El Study of the 
Hubbard model along these lines is currently underway. 
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FIG. 4. Imaginary part of the d self-energy components So, 
El, and E2, defined from the expansion Eg = Eo + Eiej + E2e| 
[see Eq. (^], and related to the local, nearest-neighbor, and 
next-nearest-neighbor self-energies through E^ — Eo -I- E2, 
E<i,j> = T,i/V2D, and E«ij» = Ai'«ij»E2/2_D [see 
Eqs. (0)]. Here Uf + Ud = 1, U/ti = 4, T/ti = 1, and D = 2. 
In the physically relevant regime, jt2/ii| ^ 1, the imaginary 
parts of El and Eq are comparable in size, indicating that im- 
portant nonlocal dynamical contributions are omitted in the 
DMFT upon neglecting t2. 
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